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ABSTRACT 

Large inhomogeneities in neutral hydrogen in the universe can be detected at redshifts 
z < 10 using the redshifted 21cm line emission. We use cosmological N-Body simula- 
tions for dark matter and a simple model for baryonic collapse to estimate the signal 
expected from structures like proto-clusters of galaxies at high redshifts. We study : (i) 
the standard CDM model, (ii) a modified CDM model with less power at small scales, 
and (iii) a A+CDM model in a universe with CIq + fl\ = 1. We show that it should 
be possible for the next generation radio telescopes to detect such structures at the 
redshift 3.34 with an integration of about 100 hours. We also discuss possible schemes 
for enhancing signal to noise ratio to detect proto-condensates at high redshifts. 
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1 INTRODUCTION 

It is generally believed that large scale structures like galax- 
ies and clusters of galaxies formed from small initial inho- 
mogeneities via gravitational collapse. One implication of 
this picture is a distinct epoch when structures like proto- 
galaxies and proto-clusters decoupled from the largely ho- 
mogeneous universe. Present observations suggest that this 
epoch is around z ~ 5 for galaxies. Observations of struc- 
tures in this stage of formation, if made, can be a very pow- 
erful constraint on the models of structure formation. Such 
observations will also improve our understanding of the pro- 
cess of 
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Su nyacv and Zcl'dovich (1072) [also 3cc Sunyacv and 

Zel'dovich (1975)] pointed out that the formation of first 
structures may be probed by observing the redshifted 21cm 
line emitted by the neutral hydrogen in these structures. 
Several searches have been made to look for such structures 
at high redshifts. In absence of a detection, these searches 
have only been able to put limits on the mass of neutral 
hydrogen present in clumped form. For a summary of these 
surveys see Wieringa, de Bruyn and Katgert (1992) and ref- 
erences cited in that paper. 

The Giant Meter-wave Radio Telescope (GMRT) 
presently being constructed in India should be able to im- 
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prove the observational situation considerably as regards the 
detection and study of proto condensates containing neu- 
tral hydrogen flSwarup 198jj ). The GMRT will be able to 
probe the redshifted 21cm line from three epochs centred 
at z = 3.34, 5.1 and 8.5. In this paper we will discuss the 
possibility of detection at the two lower redshifts. 

Subramanian and Padmanabhan (1993) have computed 
the expected flux at these redshifts for some models of struc- 
t urc formation. They used the Press-Schechter formalism 
(Press and Schechter 1975) to compute the expected num- 
ber densities of proto-clusters in the CDM and HDM mod- 



els. In a la ter paper (Kumar, Padmanabhan and Subrama- 



nian 1995) they computed line profiles assuming the proto- 



clusters to be spherically symmetric. These studies suggest 
that it should be possible to detect proto-clusters in the 
standard CDM model using the GMRT with 10 to 20 hours 
of observations. 

In this paper we will use N-Body simulations to follow 
gravitational collapse of dark matter and some simple ap- 
proximations for estimating the collapsed and neutral frac- 
tion in dense regions to construct "radio maps" with same 
specifications as the GMRT. We will study these maps and 
suggest simple methods to optimise the signal to noise ratio. 

Some authors have studied the distribution of neutral 
hydrogen at high redshifts using simulations that include 
gas dynamics, ionisation and other astrophysical processes. 
Most of these studies focus on small scale yariations in the 
distr ibution of neutral hydrogen, [e.g. see (Weinberg et al 



199c)] However, the synthesised beam for the central square 
of the GMRT includes a large comoving volume in each pixel 
(approximately 8(fa -1 Mpc) 3 ) and so the details of physical 
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processes operating at small scales can be ignored to a large 
extent. We can also ignore the differences in distribution 
of baryons and dark matter at small scales. Further, as we 
are interested in the neutral fraction at two epochs, we can 
choose to ignore the physical processes responsible for its 
evolution. This simplifies the problem to a large extent and 
we should be able to get meaningful estimates of the sig- 
nal strength without a detailed treatment of baryons and 
astrophysical processes. 

A brief discussion of the expected physical conditions at 
the epochs of interest is given below. These are the guiding 
considerations in choosing the simplifying assumptions for 
constructing the radio maps. 

z — 3.34 : Observations show th at the inter-galact 



med ium (fGM) is completely ionised (Giallongo E. et al 
this redshift. Numerical simulations suggest that 
nearly all of the neutr al hydrogen is in high density, radia- 
tively cooled objects ( Weinberg et al. 1996| ). Observations 
show that a large fraction of mass in the damped lyman- 
a absorption systems (DLAS), believed to be proge nitors of 
present day galaxies, is in form of neutral hydrogen flWolfe et| 
al. 199^). Observations also show that the spin temperature 
of gas in DLAS, the relevant quantity for the 21cm transi- 
tion, is much higher than the temperature of the background 
radiation at this epoch. It can be shown that, in such con- 
ditions, emission is the dominant mechanism and the total 
energy emitted by radiato rs in such a state doe s not depend 
on the spin temperature (Scott and Rees 1990). 

• z = 5.1 : We have very little information about the 
universe at this redshift. It is known that the intergalactic 



mediu m at redshifts z < 5 is fully ionised (Giallongo E 



et al. 1 J94 ) . It is believed that the process of ionisation is 
initiated by the first luminous structures in the universe. 
First structures like galaxies and quasars can form after z ~ 
5 — 6 in most models that satisfy other observational bounds. 
There are three possible scenarios that deserve mention. 

- One possibility is that the universe was already 
reionised by z = 5.1. In this case the neutral hydrogen 
would be confined to dense regions like proto-galaxies 
and the intervening regions will be completely, or mostly, 
ionised. 

- If the universe has been reheated by first luminous 
objects, but not reionised, then the spin temperature will 
be much higher that the temperature of the background 
radiation. A patchy reheating could lea d to fluctuations 



in the 21cm em ission at large scales (Madau, Meiksin 
and|Rees 19971). 



- If the universe is neither reheated not reionised by 
z = 5.1 and the first luminous objects form around this 
epoch, then the spin temperature will be comparable to 
the temperatur e of the background r adiation for models 
with D. b < 0.1 flScott and Rees 199C| ). 



We will focus on the second scenario (for z = 5.1) and 
carry out all our calculations assuming that the universe is 
largely neutral and the spin temperature is much larger than 
the temperature of the radiation bac kground. Fluctuations 



in spin temperature at large scales (Madau, Meiksin and 



Rees 19|97|) introduce uncertainty in the results as the region 
of interest may not have been reheated and hence the 21cm 
radiation will be much less than anticipated. 



2 GENERATING RADIO MAPS 

We will now outline the method that is used to construct 
radio maps by combining the distribution of dark matter 
obtained from N-Body simulations with approximations for 
the neutral fraction of gas. Radio observations with an in- 
terferometer give the flux of radiation coming from a given 
direction in the sky in a range of frequencies. This range 
of frequencies is subdivided into small equal intervals (chan- 
nels). Within the field of view (primary beam) we can differ- 
entiate between flux from directions separated by an angle 
equal to the resolution of the interferometer (synthesised 
beam). Therefore we can arrange the information obtained 
from such an observation in terms of a radio map for each 
frequency channel. The angular resolution, or the pixel size, 
with which the map is constructed is given by the synthe- 
sised beam of the interferometer. For the purpose of gen- 
crating "radio maps" from simulations, we will choose the 
parameters like pixel size and channel width to be s ame as 
that used by the GMRT central array ( |Swarup 1984] ). (The 
central array consists of 12 antennas of diameter 45m spread 
over a region of 1km x 1km. The telescopes are scattered 
randomly within this region.) 

Earlier estimates of expected flux from high redshift 
objects ha ve shown proto-clusters to be the mos t promis- 
ing source (Subramanian and Padmanabhan 1993) of 21cm 
radiation. Observations of high redshift objects like quasars 
suggest that the intergalactic medium is completely ionised 
at redshifts z < 4. This implies that the only source of neu- 
tral hydrogen will be dense clouds inside galaxies. These 
clouds are sufficiently dense so that absorption in a thin 
layer near the surface shields the inner regions and a large 
fraction of the gas in these remains unionised. Therefore, 
while computing the flux, we must take into account the 
fact that the neutral hydrogen we observe resides in dense 
clumps and shares the internal velocity dispersion of these. 
We take this fact into account by convolving the line profile 
with a Gaussian of width 200km/s, at the epoch of emis- 
sion. This consideration is relevant only for the window at 
z = 3.34 as in most models galaxies have not formed by 
z = 5.1. Any proto-galaxies that may exist are not expected 
to have any systematic velocity dispersion. 

A remaining uncertainty in computing the amount of 
neutral hydrogen at high redshifts is the fraction of mass 
in galaxies at the epoch of interest. We will reduce the un- 
certainty due to this factor by considering only the regions 
with S > 1 in the numerical simulations. (Changing this 
threshold to 5 > 3 reduces the peak signal by less than 
15%. Therefore the final result is not very sensitive to this 
threshold.) However, since this is a complex issue, we will 
only parameterise our ignorance with a factor f ga i which is 
the fraction of mass in galaxies that can hold neutral hy- 
drogen. We will use f ga i — 1 but the results can be rescaled 
with any value. In a more detailed calculation, which can 
be done with N-Body simulations that use a much larger 
number of particles, individual halos of mass M > 1O 1O M0 
can be directly identified and the uncertainty parametrised 
by fgai becomes irrelevant. This can not be done at present 
without reducing the physical size of the simulation box by 
a significant amount - in which case we will not be able to 
map all scales probed by the GMRT. 
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The radio map is generated from the N-Body data in 
the following manner: 

• Use the co-ordinates of the particle to fix the pixel to 
which the contribution is to be added. 

• Use the radial component peculiar velocity of the par- 
ticle to compute the total redshift. Convolve the line profile 
with a Gaussian of width 200km/s and add the contribution 
to the relevant frequency channels. 

• For generating maps at redshifts z < 4 impose a cutoff 
on density contrast 8 > 1 to ensure that we are looking at 
regions that can host galaxies. 

• Repeat this process for all particles. 

In the following subsections we will describe each of these 
components in greater detail. We begin by outlining the 
models of structure formation that are used for this study. 



2.1 N-Body Simulations 

We used N-Body simulations of three models to generate 
the "radio maps." These were chosen from the family of 
CDM models and we used the parametrised spectrum given 
by Efstathiou, Bond and White (1992). These models were 
normalised using the root mean square amplitude of fluctu- 
ations in the temperature of the CMBR (Cos mic Microwave 



Back ground Radiation) observed by COBE (Wright et al 



i99e) 



All the simulations were carried out using 128 3 particles 
in a 128 3 mesh. The physical size of the box in each case 
was 128/i _1 Mpc. Thus the mass of each particle in these 
simulations equals 2.7 x 10 QoMq. 

The following parameters were chosen for these models : 

• Model I : The standard CDM model (sCDM). (r = 
Q nr h = 0.5, Q„ r — Oo = 1 and h — 0.5.) Normalisation with 
fluctuations in the CMBR gives a(8ft- 1 Mpc,« = 0) ~ 1.2 
which is about twice the value implied by the abundance of 
rich clusters. This model is a reference model in studies of 
structure formation. 

• Model II : A flatter version of the CDM model with 
T = 0.3, Q„ r = f2o = 1 and h — 0.5. The choice of T for this 
model is independent of the cosmological parameters. This 
model has a(8h~ Mpc, z = 0) ~ 0.6 which is consistent 
with the observed amplitude of fluctuations at this scale. 
This model also predicts the correct slope of the correlation 
function at large scales. At large scales this model is a close 
approximation for some "flat" versions of CDM that can 
arise if a small but non-negligible fraction of mass in the 
universe is contributed by relativistic dark matter. (Mixed 
Dark Matter or MDM hereafter.) 

• Model III : A A+CDM model (LCDM). (r = Q nr h = 
0.3, Cl nr = 0.6, Ov = 0.4 and h = 0.5.) This model has 
cr(8/i -1 Mpc, z = 0) ~ 1 which is also consistent with the 
observed abundance of rich clusters. [The observed ampli- 
tude of fluctuatio ns at the scale of 8ft" 1 Mpc scales a pproxi- 
mately as fi~°' 6 ([White, Efstathiou and Frcnk 1992]) .] This 



model satisfies most observational constraints that are avail- 
able for cosmological models and models of stru cture forma- 
tion. (Bagla, Padmanabhan and Narlikar 1996) (These con- 
straints rule out large regions of the parameter space and 
only a small region survives. Therefore model III is a good 





Figure 1. This figure shows the power spectra for the three mod- 
els being discussed here. The standard CDM model is shown us- 
ing a thick line, model II (MDM) as a dashed line and model III 
(LCDM) as a dot-dashed line in these panels. We have shown the 
linearly evolved power spectrum at the present epoch [top panel] 
and also at the redshift z = 3.34 [lower panel]. The different rate 
of growth for model III (LCDM) in the linear regime leads to 
an increase in the ratio of power in this model to the power in 
model I (sCDM) or model II (MDM). 



example from the class of allowed models.) This model was 
chosen for two reasons : 

- Growth of perturbations in a A model slows down 
at late epochs. In other words, such a model has more 
power at early times in comparison with an fi nr = 1 
universe which has same level of clustering at z = 0. 
Thus we have reason to expect higher signal for models 
in a universe with nonzero A. 

- The comoving volume enclosed in a given solid an- 
gle at high redshifts is higher for a universe with nonzero 
A. This implies more emitters and hence a higher signal. 
In a model with nonzero cosmological constant the lu- 
minosity distance is also larger in comparison with the 
corresponding distance in the Einstein-de Sitter model 
for any given redshift. However, the increase in comoving 
volume mentioned above compensates for the increase in 
luminosity distance in such models. 

We have shown the power spectra for these models in fig- 
ure 1. To demonstrate the effect of different rate of growth 
in linear regime we have also shown the same spectra at 
redshift 3.34. 

We would like to emphasise the following points regard- 
ing the choice of models: It is possible to choose many vari- 
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ants of the above models and repeat the analysis given in 
this paper. One can choose to normalise the power spectra 
in a different manner, e.g. by fixing the amplitude of per- 
turbations at 8h~ Mpc rather than at very large scales us- 
ing COBE observations. Other parameters like the Hubble's 
constant, the cosmological constant, the density parameter 
and the primordial spectral index can also be varied from 
the values used here. We have consciously restricted our- 
selves to a small set of example models for the purpose of 
illustration. A more detailed exercise will be warranted if 
the GMRT succeeds in detecting neutral hydrogen at these 
redshifts and a more meaningful comparison with models 
can be carried out using real data. 



2.2 Evolution of Neutral Fraction 

Generating artificial radio maps requires, in principle, a de- 
tailed knowledge of the distribution of baryons and the neu- 
tral fraction etc. However, as mentioned earlier in this paper, 
the comoving volume enclosed in each pixel is very large and 
therefore we can get reasonable estimates by assuming the 
distribution of baryons and dark matter to be the same. In 
the following discussion we will demonstrate, using a simple 
model, that the assumption of a constant neutral fraction at 
these scales is not very far from the truth. We will also show 
that our choice of /jv = 0.5 for z = 3.34 is not unrepresenta- 
tive. We will show that the neutral fraction of gas in galaxies 
is largely independent of the depth of the potential well. In 
the following discussion we estimate the neutral fraction by 
using the model of star formation in galaxies of Kauffmann, 
White and Guiderdoni (1993; hereafter KWG). Here we will 
briefly summarise the relevant features of their model and 
use these to estimate the evolution of neutral fraction. 

In this model, dark matter halos are assumed to be 
truncated singular isothermal spheres and it is assumed that 
the temperature T of the gas is given in terms of the circular 
velocity, V c , as, T = 35.9 (V c /km/s) 2 K. The virial radius r v 
is defined to be the radius within which the mean overdensity 
is 200 (i.e., r v = 0.1H _1 (1 + z)~ 3/2 V c ). The radius where 
the cooling time of the gas is equal to the age of the universe 
is defined as the cooling radius, r coo i. 

Suppose that the fraction of the critical density that 
is in baryons is fi& and f g is the fraction of the baryons 
in the form of gas. The amount of cold gas inside the halo 
at time t is given by the amount of gas with cooling time 
tcooi < t. When r coo i 2> r v , cooling is very rapid and the 
rate of increase of cold gas in the halo is governe d by the 
accretion rate of the halo, (White and Frenk 1991) 



Minf(V c , z) = 0.15 f g Q b V? G- 1 



(1) 



In the other limit, r coo i <C r v , the rate of inflow of cold gas 
can be written as, 



M cool {V c ,z) = ATYp g (r cool )r cool - 



dt 



(2) 



where p g (r) is the gas density at radius r. 

In the model of KWG, the rate at which cold gas settles 
inside the halo is given by min(M in /, M coo i). For the cooling 
of the gas, we use the cooling function of Fall and Rees 
(1985) for a primordial gas, since the metallicity is any way 
small for gas at high redshift. 



However, with the onset of star formation, the super- 
novae will begin to heat the gas. Following KWG, we assume 
that the number of supernovae per solar mass of stars formed 
is t]sn — 4 x 10 -3 Mg 1 . If each of the supernovae have a 
kinetic energy of the ejecta of order 10 51 erg, and a fraction 
e of this energy is used to heat the cold gas, then the rate of 
loss of cold gas to the hot phase of the int erstellar medium 
is, ( |Kauffman, White and Guideroni 1993 ) 



Mj 



4 M«rjsNE S N 

cheat — £ g y 2 



(3) 



Here, M» is the star formation rate, which is given in this 
model as, 

M coU 



M, 



tdy 



(4) 



Here tdyn is the dynamical time. The authors defined the dy- 
namical time as (where A ~ 0.05 is the initial dimensionless 
spin parameter of the gas), 



tdyn 



V gal ~ (GM gal y/i 



(5) 



These equations govern the evolution of the amount of 
cold gas, or, in other words, the neutral fraction fni, of gas 
in a halo of circular velocity V c at a given redshift z, for an 
assumed value of the onset of inflow of gas z;„, and given 
the values of a and e. Gas is assumed to be fully ionised 
at Zi n . KWG estimated the values of a and e from evolving 
the stellar population according to the above model and 
comparing the mean luminosity and the cold gas content 
of halos of Vc = 220 km/s with the observed values for 
the Milky Way [their Table 1]. For example, for fio = 1, 
fib = 0.1, a = 0.1, and e = 0.08. The values of a and e are 
smaller for smaller values of fib. However, we have found 
that the neutral fraction is not very sensitive to the values 
of a and e, but rather to z; n . 

We plot in figure 2 the neutral fraction /hi at z = 3.34 
for Zi n = 4, 5, 6, 9, as functions of the rotational velocity V c . 
It is seen that the neutral fraction is not sensitive to the 
value of V c , for V c > 150 km/ s. We will, therefore, use a 
constant value of /„ for all potential wells. We also show 
the values of /hi at z = 5.1 for Zi n = 6, 9 as functions of 
V c in figure 2. The neutral fraction f n does not depend very 
strongly on the value of / 9 fii,. In the following discussion we 
will use /„ = 0.5. 

2.3 Redshift Space Projection 

Simulations output the data in form of positions and veloc- 
ities of particles whereas radio observations detect the flux 
as a function of angular position and frequency. In this sub- 
section we will outline the method used for computing the 
observed central frequency around which the 21cm radia- 
tion from hydrogen atoms represented by a given N-Body 
particle is observed. 

The peculiar velocity of an N-Body particle is given by 



dx 

da 



H a 



H(a) dx 



(6) 



Ho da 

We only need the radial component of velocity for comput- 
ing the redshift. As we are dealing with particles at very 
high redshifts we can align one axis of the simulation box 
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a combination of the Hubble recession and the peculiar ve- 
locity. 



1 + Ztot = (1 + Z hub )(l + Z pec ) 



(8) 



The average z^ub for the simulation box is 1/a—l. Redshift 
due to peculiar velocity is given by z pec = v p /c where c is 
the speed of light. 

2.4 Flux 

If the spin temperature of the atoms emitting 21cm radiation 
is much greater than the temperature of the CMBR then the 
spin temperature drops out of the expression for the emitted 
energy. In such a case the energy dE emitted by a set of 
hydrogen atoms in an interval <ft e is given by 

dE — Rate of Transition x Energy carried by a photon 
x Number of hydrogen atoms x dt e 

(9) 



3 A U M Hj M 
-A2I X hv e X -dt e 

4 m v 



where m p is the mass of a proton and Mh i is the total mass 
in neutral hydrogen. We will use the mass in neutral hydro- 
gen contributed by one particle in the N-Body simulation to 
compute total energy emitted by each "N-Body particle". 
The mass Mh, is given by 



Ml 



Mpartfbfn 



= 8.1 X 10 Mf 



M ^ b f 

mpart J n 

"0 

( SH \ (fn\ 
VO.06/ \0.5J 



(10) 



here /(, is the fraction of baryons, fif, is the contribution 
of baryons to the density parameter and f n is the neutral 
fraction. We have chosen fit = 0.06 as this value compares 
well with the observed a bundance of light elements and pr i- 
mordial nucleosynthesis ( Copi, Schramm and Turner 1995 ) . 
Using this, we can estimate the flux received by an observer 
from an "N-Body particle". For Q nr = f2o = 1 models the 
flux contributed by one particle at redshift z — 3.34 is given 
by 



S v = 1.1/xJy 



M 



8.1 x 10 9 M Q I \l75kHz 



Hi 



(11) 



The frequency width used here corresponds to a velocity 
dispersion of 200km/s. (The corresponding number for the 
model with Q A = 0.4 is 0.77fiJy.) 



Figure 2. The upper panel shows neutral fraction fu j at 2 = 
3.34 against the circular velocity V c of the halo for «j n = 4, 5, 6, 9. 
We have assumed a = 1, e = 0.08, fgflt, = 0.05, h = 0.5 and 
fi = 1. Here, f g is the fraction of baryonic mass in the form of 
gas. The curves do not strongly depend on the value of f g . The 
lower panel similarly plots the neutral fraction at 2 = 5.1 for 

Zin — 6, 9. 

along the line of sight and use that component of velocity 
for computing the recession velocity. The ratio of H/Ho is 
given by 

for the models of interest. The effective redshift results from 



3 RESULTS 

In this section we shall outline the results of analysis of the 
radio maps generated from N-Body simulations. We begin 
with a pictorial preview of the radio maps. Figure 3 shows 
a sample radio map for each of the three models at red- 
shift z = 3.34. The panels of this figure show one frequency 
channel (chosen to be 125kHz; which - in velocity units - 
corresponds to about 115km/s.) The contours in these ra- 
dio maps correspond to 15, 30 and 60/^Jy. The pixel size is 
3.2 arc minutes and it corresponds to a comoving scale of 
2.7/!~ 1 Mpc for models I and II. (3.5/i _1 Mpc for model HI.) 
It is clear from these panels that models I (sCDM) and III 
(LCDM) have comparable signal whereas model II (MDM), 
as it has less power at smaller scales, has somewhat lower 
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Figure 3. This figure shows a sample radio map for model I 
(sCDM) at redshift z = 3.34. This map is shown for one frequency 
channel of width 125kHz [corresponding to a velocity width of 
115km/s]. Each pixel corresponds to an angular resolution of 3.2 
arc minutes and a physical scale of about 3/i _1 Mpc. The contour 
levels correspond to 15, 30 and 60/^Jy. 




Figure 3. Continued. This figure shows a sample map for model 
II (MDM). 



signal. The profile of peaks in the sCDM model is steeper 
than that in the LCDM model. 

Figure 4 shows similar maps for z = 5.1 \y = 233MHz]. 
The width of one channel in this case corresponds to a ve- 
locity width of 161km/s. Angular size of each pixel is 4.5 
arc minutes and this corresponds to a comoving scale of 
4.7ft~ 1 Mpc for models I (sCDM) and II(MDM). (5.4/i _1 Mpc 




Figure 3. Continued. This figure shows a sample map for model 
III (LCDM). 



for model III (LCDM).) These maps clearly show that there 
are no small scale structures in these, in comparison with 
the maps at z = 3.34. One reason for the gentler variation 
in signal from one pixel to another arises from the fact that, 
having assumed that the IGM is not ionised, we have not dis- 
carded matter in under dense regions while computing the 
signal. The signal is much higher as the comoving volume 
enclosed in each pixel/channel is larger than the correspond- 
ing volume at redshift 3.34. 

These qualitative features can also be seen in figure 5 
that shows a sample spectrum from the simulations of 
model III (LCDM) at these two redshifts. The spectra are 
shown as a function of channel number and the signal is 
shown in fijy. The typical width of high peaks seen in 
these graphs are much larger than the velocity dispersion 
of 200km/s ascribed to individual particles. 

In order to make quantitative estimates of the possibil- 
ity of detecting neutral hydrogen at high redshifts we com- 
puted the amplitude of three highest peaks in the radio maps 
for the three models. Table 1 lists the amplitudes and the 
full width at half maximum for these models for the three 
highest peaks in each case at redshift 3.34. The last column 
lists the number of pixels enclosed within the contours of 
half maximum. The contours were drawn for signal aver- 
aged over frequency channels. The number of channels used 
for averaging was taken to be the FWHM for the relevant 
peak. This number indicates the typical angular size of a 
flux peak and the signal can be smoothed at this scale to 
enhance signal to noise ratio. However, the shapes of these 
peaks are somewhat arbitrary and it is difficult to suggest a 
generic smoothing function for improving chances of detec- 
tion. Figure 6 shows contours around one such peak, showing 
the level at which the signal drops to .75, .5 and .25 its max- 
imum value. This figure demonstrates the point mentioned 
above regarding shapes of these contour. 

Even though we cannot suggest the most optimum 
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Figure 4. This figure shows a sample radio map for model I 
(sCDM) at redshift z = 5.1. This map is shown for one frequency 
channel of width 125kHz (corresponding to a velocity width of 
160km/s). Each pixel corresponds to an angular resolution of 4.5 
arc minutes and a physical scale of about 5h~ x Mpc. The contour 
levels correspond to 40, 80, 120 and 200/^Jy. 




Figure 4. Continued. This figure shows a sample map for the 
MDM model (model II). 



smoothing function for improving the signal to noise ra- 
tio, we have confirmed that the signal to noise ratio does 
improve for smoothing with a square top hat window of 

1/2 

size w n^ M where uhm is the number of pixels enclosed 
within the contour of Smax/2. We assumed that the noise 
scales as n~ 1,/2 for small n. The signal to noise ratio for mod- 




Figure 4. Continued. This figure shows a sample map for the 
LCDM model (model III). 



Model 


Amplitude (fiJy) 


FWHM (MHz) 


riHM 


I 


137.2 


1.5 


3 




101.5 


1.0 


11 




99.8 


0.875 


7 


II 


87.1 


1.0 


6 




97.7 


0.75 


8 




93.5 


0.875 


12 


III 


124.6 


0.875 


11 




111.3 


0.875 


13 




103.8 


1.125 


8 



Table 1. This table lists details of three highest peaks in the 
simulated radio maps for the three models of structure forma- 
tion being considered here for the redshift window at z = 3.34 
(327MHz). We have listed the amplitude of peaks in micro Jan- 
skys, full width at half maximum in MHz and the number of pixels 
enclosed in the contours of half maximum. It is clear that mod- 
els I (sCDM) and III (LCDM) have a higher signal as compared 
to model II (MDM). 



two in this manner. The standard CDM model [model I] 
has more small scale power and hence the peaks are much 
sharper and therefore the gain in signal to noise ratio by 
smoothing is somewhat limited in this case. 

In order to estimate the integration time required for 
detection and imaging of these peaks we need to know the 
root mean square amplitude of noise expected for the GMRT 
receivers. The expected noise for the central array is 



rms noise = 



4, r f T a \ flMHz\ 1/2 /1007ws\ 1/2 
WOfiJy ( 



HOif 

T, 



250K 



Av 

nMHz \ 1/2 
\ Av ) 



The system temperature for 327MHz window is 11 OK and 
the c orresponding number for 233MHz is 250K. (Swarup 



els II ( tylDM) and HI (LCDM) can be improved by a factor 198j ) 
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Figure 5. This figure shows a sample spectrum for each of the 
redshifts. The spectra shown here are taken from the simulated 
radio maps of the LCDM model (model III). The top panel cor- 
responds to z = 3.34 and the lower panel to z = 5.1. The sig- 
nal in /ijy is shown as a function of frequency channel [channel 
width=125KHz] . This figure shows that the signal at z = 5.1 
has few small scale variations in comparison with the signal at 
z = 3.34. 

It is clear from table 1 that the peak signal expected 
in models I (sCDM) and III (LCDM) is about three times 
larger than the rms noise expected at redshift 3.34. Further, 
the width of the peaks in frequency space is comparable to 
1MHz. The signal to noise level can be enhanced further by 
smoothing over nearby pixels. The number of pixels enclosed 
in the contour of S ma x/2 for model III (LCDM) varies be- 
tween 8 and 13 for the three higher peaks listed in table 1. 
Therefore an enhancement by a factor two or more can be 
obtained by smoothing over nearby pixels. This implies that 
a 3<r detection of such objects should be possible with an 
integration time of 50 — lOOhrs. The main uncertainty in 
this result comes from our lack of knowledge about the frac- 
tion of mass in galaxies. We have, of course, tried to reduce 
the uncertainty by using only regions with density contrast 
greater than unity for computing the signal. As mentioned 
before, we find that changing the density threshold to 5 C = 3 
does not change the amplitude of the highest peaks by more 
than 15%. However, as the estimated signal and noise have 
a similar amplitude, even a small factor may make all the 
difference. 

Detection of neutral hydrogen at high redshift for 
model I (sCDM) should also require similar integration time. 




Figure 6. This figure shows the contours for / * Smax around 
a peak for / = 0.75, .5 and 0.25. These contours is for one of 
the highest peaks in simulated maps for model III (LCDM). The 
peak is located at (5,5) in this picture and is marked by a star. 
It is clear from this figure that the signal can be averaged over 
neighbouring pixels to improve signal to noise ratio. However the 
shape of of this contour suggests that the amount of gain from 
such an optimisation will be limited. 



Here the frequency spread of the highest peak is very large 
and compensates, to some extent, the compactness of the 
peak in angular co-ordinates. On the other hand detection 
of neutral hydrogen in model II (MDM) is a more difficult 
proposition as the peak signal is small compared to that 
for the other two models. The integration time required for 
model II (MDM) is at least twice as large as that for models I 
(sCDM) and III (LCDM). 

Table 2 lists the amplitude of three highest peaks for the 
three models we are using here for the window at redshift 
5.1. We have also listed the full width at half maximum for 
each of these peaks. These peaks have a very shallow profile 
and the number of pixels, enclosed within the region where 
the signal is greater than half of the maximum value, can be 
very large. Therefore we can smooth the signal over a few 
neighbouring pixels and improve the signal to noise ratio 
by a significant amount. Results for different models in this 
case are: 

• Model III (LCDM) predicts a signal that is comparable 
to the rms noise. Smoothing over nearby pixels and integra- 
tion for 50 — lOOhrs should be sufficient for a 3a detection. 

• Model I (sCDM) also predicts signal at the same level 
as model III (LCDM). 

• Model II (MDM) has very low signal as compared to 
models I (sCDM) and III (LCDM) and will require an inte- 
gration for 100 — 200 hours for a 3<r detection. 

These numbers, however, must be considered in light 
of the possibility of a patchines s in reheating of the IGM. 
( Madau, Meiksin and Rees 1997 ) 
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IVlOtlUl 


Amplitude (^fijy^j 


r vv inivi i ivinz t 


1 


199.0 


0.75 




181.4 


1.0 




1 fV^ 9 




TT 
11 


loo.O 


U. f 




115.2 


0.875 




103.8 


0.75 


III 


221.7 


0.75 




195.6 


0.875 




153.9 


0.875 



Table 2. This table lists details of three highest peaks in the sim- 
ulated radio maps for the three models of structure formation be- 
ing considered here for the redshift window at z = 5.1 (233MHz). 
We have listed the amplitude of peaks in micro Janskys and full 
width at half maximum in MHz. It is clear that models I (sCDM) 
and III (LCDM) have a higher signal as compared to model II 
(MDM). 



4 DISCUSSION 

Results of numerical simulations suggest that it could be 
possible to detect neutral hydrogen at high redshifts, in 
proto-clusters. Detection of these objects with the GMRT 
will require integration over 50 — 100 hours. The signal to 
noise ratio for z — 3.34 is better than that expected for 
z = 5.1 if we choose the factor f ga i to be unity. This is the 
most uncertain number in our calculation. If this number 
is much smaller than unity then it may be difficult to im- 
age proto-clusters at this redshift. It may still be possible to 
"detect" neutral hydrogen statistically by doing a (AT)/T 
type of an experiment. In these one looks for excess corre- 
lations than are expected from noise in the instrument. For 
example, the emission from a region of size 10 6 (/i _1 Mpc) 3 
- typical size of the volume probed in one field of view - is 
around 370mJy. [This number is obtained by using the av- 
erage flux in the simulated radio maps.] This will appear as 
excess noise if the individual structures can not be imaged. 
Integration over more than one field of view can also be used 
to enhance the chances of de tection as one may be able to 



pick up a rare density peak ( subramanian and Padmanab- 
han 1993| ). 

As mentioned above, the biggest uncertainty in our re- 
sults for z — 3.34 is introduced by our lack of knowledge 
about f g ai ■ This uncertainty can only be removed with large 
simulations that can resolve masses less than 1O 1O M0 and 
still cover a volume that is comparable with that covered 
by one field of view for the GMRT. The other option is to 
carry out a series of N-Body experiments at different scales 
for each of the models and deduce f ga i for each model from 
the ensemble of simulations. 

A c omparison with earlier estimates for the expected 
signal (Kumar, Padmanabhan and Subramanian 199E) 



shows that the highest signal obtained in simulations is 
smaller than that expected from a proto-cluster of mass 
1O 15 M0. This could result due to three reasons 

• We do not have a sufficiently high peak in our realisa- 
tion of the density field. 

• The angular extent of the proto-clusters is considerably 
larger than one pixel. 



• The velocity dispersion of these clusters is larger than 
expected. 

We feel that all of these factors have contributed to make 
our estimates smaller than the earlier estimates by a factor 
2-5. 

The possibility of detection of neutral hydrogen at red- 
shift 5.1 may be easier if the universe has been reheated but 
not fully reionised. Another good feature of 21cm emission 
from this epoch is that it only has smooth large scale vari- 
ations in signal and hence it should be possible to enhance 
the signal to noise ratio by at least a factor 2 — 3 by smooth- 
ing over nearby pixels. However, the results for this epoch 
are somewhat uncertain due to our poor knowledge of the 
physical conditions at z > 5. 
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